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Abstract 

In a metastable de Sitter space any object has a finite life expectancy beyond which it 
undergoes vacuum decay. However, by spreading into different parts of the universe which 
will fall out of causal contact of each other in future, a civilization can increase its collective 
life expectancy, defined as the average time after which the last settlement disappears due to 
vacuum decay. We study in detail the collective life expectancy of two comoving objects in de 
Sitter space as a function of the initial separation, the horizon radius and the vacuum decay 
rate. We find that even with a modest initial separation, the collective life expectancy can 
reach a value close to the maximum possible value of 1.5 times that of the individual object 
if the decay rate is less than 1% of the expansion rate. Our analysis can be generalized to 
any number of objects, general trajectories not necessarily at rest in the comoving coordinates 
and general FRW space-time. As part of our analysis we find that in the current state of the 
universe dominated by matter and cosmological constant, the vacuum decay rate is increasing 
as a function of time due to accelerated expansion of the volume of the past light cone. Present 
decay rate is about 3.7 times larger than the average decay rate in the past and the final decay 
rate in the cosmological constant dominated epoch will be about 56 times larger than the 
average decay rate in the past. This considerably weakens the lower bound on the half-life of 
our universe based on its current age. 
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1 Introduction 

The possibility that we may be living in a metastable vacuum has been explored for more that 
hfty years [1-6] . Discovery of the accelerated expansion of the universe [7,8] and subsequent 
developments in string theory leading to the construction of de Sitter vacua [9-11] suggest that 
the vacuum we are living in at present is indeed metastable. Unfortunately our understanding 
of string theory has not reached a stage where we can make a definite prediction about the 
decay rate of our vacuum. The only information we have about this is from the indirect 
observation that our universe is about 1.38 x 10^° years old. Therefore, assuming that we have 
not been extremely lucky we can conclude that our inverse decay rate^ is at least of the same 

^For exponential decay the inverse decay rate differs from half-life by a factor of In 2. In order to simplify 
terminology, we shall from now on use only inverse decay rate and life expectancy - to be defined later - as 
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order. ^ 

Typically the decay of a metastable vacuum proceeds via bubble nucleation [1-6] (see [14] 
for a recent survey). In a small region of space-time the universe makes transition to a more 
stable vacuum, and this bubble of stable vacuum^ then expands at a speed that asymptotically 
approaches the speed of light, converting the rest of the region it encounters also to this stable 
phase. Due to this rapid expansion rate it is impossible to observe the expanding bubble before 
encountering it - it reaches us when we see it. However, due to the existence of the future 
horizon in the de Sitter space, even a bubble expanding at the speed of light cannot hll the 
whole space at future inhnity. Indeed, it has been known for quite some time that in de Sitter 
space if the expansion rate of the universe exceeds the decay rate due to phase transition then 
even collectively the bubbles of stable vacuum cannot £11 the whole space [15] and there will 
always be regions which will continue to exist in the metastable vacuum. Nevertheless, any 
single observer in the metastable vacuum will sooner or later encounter an expanding bubble 
of stable vacuum, and the probability of this decay per unit time determines the inverse decay 
rate of the observer in the metastable vacuum. 

This suggests that while any single observer will always have a limited average life span 
determined by the microscopic physics, a civilization could collectively increase its longevity 
by spreading out and establishing different civilizations in different parts of the universe [16]. 
If the bubble of stable vacuum hits the civilization - henceforth refered to as object - in the 
initial stages of spreading out then it does not help since the same bubble will most likely 
destroy all the objects. However, with time the different objects will go outside each other’s 
horizon and a single bubble of stable vacuum will not be able to destroy all of them. This will 
clearly increase the life expectancy of the objects collectively - defined as the average value of 
the time at which the last surviving object undergoes vacuum decay - although there will be 
no way of telling a priori which one will survive the longest. A simple calculation shows that 
if we could begin with 2 objects already far outside each other’s horizon so that their decay 

measures of longevity. 

^We shall in fact see in §5.3 that the actual lower bound for the current inverse decay rate is weaker by a 
factor of 3.7, making it comparable to the time over which the earth will be destroyed due to the increase in 
the size of the sun. Allowing for the possibility that we could have been extremely lucky reduces the lower 
bound on the inverse decay rate by about a factor of 10 [12, 13]. 

^We shall refer to the more stable vacuum as the stable vacuum, even if this vacuum in turn could decay to 
other vacua of lower energy density. In any case since this vacuum will have negative cosmological constant, the 
space-time inside the bubble will undergo a gravitational crunch [6] . We shall ignore the possibility of decay to 
Minkowski vacua or other de Sitter vacua of lower cosmological constant since the associated decay rates are 
very small due to smallness of the cosmological constant of our vacuum. 
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Figure 1: The figure showing the ‘gain’ in the life expectancy for two objects compared to that 
of one object as a function of T for r = .0003, .001, .003, .01, .03, .1 and .3. 

probabilities can be taken to be independent, then the life expectancy of the combined system 
increases by a factor of 3/2 compared to the life expectancy of a single isolated object. In the 
case of n copies the life expectancy increases by a factor given by the n-th harmonic number. 
However, in actual practice we cannot begin with copies of the object already outside each 
other’s horizon. As a result the increase in the life expectancy is expected to be lower. 

The goal of this paper will be to develop a systematic procedure for computing the increase 
in the life expectancy of the object as a result of making multiple copies of itself. For two 
objects we obtain explicit expression for the life expectancy in terms of three parameters: the 
Hubble constant H of the de Sitter space-time determined by the cosmological constant, the 
vacuum decay rate or equivalently the life expectancy T of a single isolated object and the 
initial separation r between the two objects. In fact due to dimensional reasons the result 
depends only on the combination HT and Hr, so we work by setting H = 1. In Fig. 1 we 
have shown the result for the ratio of the life expectancy of two objects to that of a single 
object - called the ‘gain’ - as a function of T for different choices of r. From this we see that 
even for a modest value of r = 3 x 10“^ the gain in the life expectancy reaches close to the 
maximum possible value of 1.5 if T is larger that 100 times the horizon size of the de Sitter 
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Figure 2: Growth of the relative decay rate R - dehned as the ratio of the decay rate to 
its asymptotic value - of a single object in FRW space-time as a function of the scale factor 
a. Value of i? at a = 1 represents the decay rate today relative to what it would be in the 
cosmological constant dominated epoch. 

space, i.e. the decay rate is less than 1% of the expansion rate. T = 100 corresponds to about 
1.7 X 10^^ years, r = 3 x 10“^ corresponds to a physical distance of the order of 5 x 10® light 
years and is of the order of the minimal distance needed to escape the local gravitationally 
bound system of galaxies. If T = 10 - i.e. of order 1.7 x 10^^ years - the gain is about 20% for 
r = 3 X 10“^. These time scales are shorter than the time scale by which all the stars in the 
galaxy will die. Therefore, if T lies between 10^^ years and the life span of the last star in the 
local group of galaxies which will be gravitationally bound and will remain inside each other’s 
horizon, then we gain a factor of 1.2 - 1.5 in life expectancy even by making one additional 
copy of the object at a distance larger than about 10^ light years from us. On the other hand 
if T is larger than the life span of the last star in the galaxy then our priority should be to 
plan how to survive the death of the galaxy rather than vacuum decay. Some discussion on 
this can be found in [17]. 

Even though most of our analysis focusses on the case of a pair of objects in de Sitter space 
at hxed comoving coordinates, our method is quite general and can be applied to arbitrary 
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number of objects in a general FRW metric moving along general trajectories. We discuss these 
generalizations in §5. In particular considering the case of a single object in an FRW metric 
dominated by matter and cosmological constant, as is the case with the current state of our 
universe, we hnd that the vacuum decay rate increases as a function of time due to accelerated 
expansion of the volume of the past light cone. This has been shown in Fig. 2. This rate 
approaches a constant value as the universe enters the cosmological constant dominated era, 
but we hnd for example that this asymptotic decay rate is about 15 times larger than the 
decay rate today, which in turn is about 3.7 times larger than the average decay rate in the 
past. Now given that the universe has survived for about 1.38 x 10^° years, we can put a lower 
bound of this order on the inverse of the average decay rate in the past.^ This translates to a 
lower bound of order 3.7 x 10® years on the inverse decay rate today and 2.5 x 10® years on 
the asymptotic inverse decay rate. 

The rest of the paper is organised as follows. In §2 we describe the case of the decay of n 
objects assuming that their decay probabilities are independent of each other, and show that 
the life expectancy of the combined system goes up by a factor equal to the n-th harmonic 
number. In §3 we carry out the complete analysis for two observers in 1+1 dimensional de 
Sitter space. The hnal result for the life expectancy of the combined system can be found in 
(3.30). This is generalized to the case of two observers in 3+1 dimensional de Sitter space- 
time in §4. Eq.(4.17) together with (4.16) and (4.11) gives the probability that at least one 
of the two objects survives till time t, which can then be used to compute the life expectancy 
of the combined system using (4.18). In §5 we discuss various generalizations including the 
case of multiple observers, general trajectories and general FRW type metric. We conclude 
in §6 with a discussion of how in future we could improve our knowledge of possible values 
of the parameters r and T which enter our calculation. In appendix A we compute the time 
dependence of the decay rate for a general equation of state of the form p = wp. 

2 Independent decay 

Let us suppose that we have two independent objects, each with a decay rate of c per unit 
time. We shall label them as Ci and C 2 - If we begin with the assumption that both objects 

^One must keep in mind that this is not a strict bound since we could have survived till today by just being 
lucky. 
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exist at time t = 0 then the probability that the hrst object exists after time t is 

Pi{t) = e-^K ( 2 . 1 ) 

Therefore, the probability that it decays between time t and t + 5t is —Pi{t)5t where Pi denotes 
the derivative of Pi with respect to t, and its life expectancy, is 

nOQ poo 

P = — / tPi{t)dt = / Pi{t)dt = c~^. (2.2) 

Jo Jo 

Independently of this the probability that the second object exists after time t is also given by 
exp[—ct] and it has the same life expectancy. 

Now let us compute the life expectancy of both objects combined, defined as the average of 
the larger of the actual life time of Ci and 6 * 2 . To compute this note that since the two objects 
are independent, the probability that both will decay by time t is given by (1 — Pi(t))(l — 
P 2 (t)) = (1 — Pi(t))^. Therefore, the probability that the last one to survive decays between t 
and t + 5t is -^{1 — Pi{t)Y5t. This gives the life expectancy of the combined system to be 

/ oo d S 

t^(l-A(t))^dt = -c-^ (2.3) 

Therefore, we see that by taking two independent objects we can increase the life expectancy 
by a factor of 3/2. A similar argument shows that for n independent objects the life expectancy 
will be 

3 Vacuum decay in 1+1 dimensional de Sitter space 

Consider 1+1 dimensional de Sitter space 

ds^ = —dt^ + e^^dx^ . (3.1) 

Note that we have set the Hubble constant of the de Sitter space and the speed of light to 
unity so that all other time / lengths appearing in the analysis are to be interpreted as their 
values in units of the inverse Hubble constant. We introduce the conformal time t via 

T=-e-* (3.2) 
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Figure 3: A comoving object in de Sitter space and its past light cones at conformal times t 
and T + 5t. 


in terms of which the metric takes the form 

= T“^(—dT^ + dx^). (3.3) 

At f = 0 we have t = — 1 and comoving distances coincide with the physical distances. 

We shall use this space-time as a toy model for studying the kinematics of vacuum decay. 
We shall assume that in this space-time there is a certain probability per unit time per unit 
volume of producing a bubble of stable vacuum, which then expands at the speed of light 
causing decay of the metastable vacuum. We shall not explore how such a bubble is produced; 
instead our goal will be to study its effect on the life expectancy of the objects living in this 
space. In §4 we shall generalize this analysis to 3-1-1 dimensional de Sitter space. 


3.1 Isolated comoving object 


Consider a single object in de Sitter space at rest in the comoving coordinate x (say at x = 0), 
shown by the vertical dashed line in Fig. 3. We start at t = 0 (t = —1) and are interested in 
calculating the probability that it survives at least till conformal time t. If we denote this by 
Fo(t^) then the probability that it will decay between t and is —Pq{x)5x where ' denotes 

derivative with respect to t. On the other hand this probability is also given by the product 
of Fo(t^) and the probability that a vacuum bubble is produced somewhere in the past light 
cone of the object between t and t-|- as shown in Fig.3. The volume of the past light cone 
of this interval can be easily calculated to be 


2 6x 


da 

8 


2 
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5t. 


— OO 


(3.4) 
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Figure 4: Two comoving objects in de Sitter space separated by physical distance r at t = — 1. 


Therefore, if K is the probability of producing the bubble per unit space-time volume then the 
probability of producing a bubble in the past light cone of the object between t and t -|- is 
given by —2K5t:/x. The previous argument then leads to the equation 

F'(t) = 2Kx-^Po{x ). (3.5) 

This equation, together with the boundary condition Pq{x = —1) = 1, can be integrated to 
give 

\nPo{x) =2K\n{-x). (3.6) 

In terms of physical time t we have^ 


Po(t) = e-^^*. (3.7) 

From this we can calculate the life expectancy, dehned as the integral of t weighted by the 
probability that the object undergoes vacuum decay between t and t -|- dt. Since the latter is 
given by —Po{t)dt, we have the life expectancy 

/ oo poo 

t Po{t) dt = j Poit) dt = (3.8) 

where in the second step we have used integration by parts. We shall express our final results 
in terms of T instead of K. 

3.2 A pair of comoving objects 

Next we shall consider two comoving objects Ci and C 2 in de Sitter space separated by physical 
distance r at f = 0 or equivalently t = —1. We shall take r < 1, i.e. assume that the two 

®Note that by an abuse of notation we have used the same symbol Pq to denote the probability as a function 
of t although the functional form changes. We shall continue to follow this convention, distinguishing the 
function by its argument (t or t). Derivatives with respect to t and t will be distinguished by using Pq to 
denote T-derivative of Pq and Pq to denote t-derivative of Pq. 
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Figure 5: The past light-come of 6*2 at t = —1 and the past light cone of Ci between t and 
T -|- for T < r — 1. 


objects are within each other’s horizon at the time they are created. We denote by T’j(T) the 
probability that Ci survives at least till conformal time t for i = 1,2 and by Pi 2 (t^iW 2 ) the 
joint probability that Ci survives at least till conformal time Ti and C 2 survives at least till 
conformal time T 2 . The boundary condition will be set by assuming that both objects exist at 
T = — 1 , so that we have 


Pi(-l) = 1, P 2 (-l) = 1, t:2) = A(t: 2 ), -1) = ■ (3.9) 


First we shall calculate Pi(t) and P 2 (t^)- They must be identical by symmetry, so let us 
focus on Pi(t). The calculation is similar to that for Po(t^) above for a single isolated object, 
except that the existence of (^2 at t = —1 guarantees that no vacuum decay bubble is produced 
in the past light-come of (72 at t = — 1 , and hence while computing the volume of the past 
light cone of the Ci between t and t -|- 5t, we have to exclude the region inside the past light 
cone of (72 at T = — 1. This has been shown in Fig. 5. This volume is given by 


(5t 



da 

7-00 



2 

r + 1-'X 


for T < r — 1. 


(3.10) 


However, for t > r — 1 the past light cone of Ci between t and t -|- does not intersect the 
past light cone of (72 at t = —1 (see Fig. 6 ), and we get the volume to be 


2(5t 



—2 - for T > r — 1. 

T 


(3.11) 
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Figure 6: The past light-come of C 2 at t = — 1 and the past light cone of Ci between t and 
T -|- for T > r — 1. 


This leads to the following differential equation for -Pi(t^): 


1 dPi _] -K [- 2 /t - 2/(r + 1 - t)] for t < r - 1, 

Pi{x) dx \2Klx for T > r — 1. 

Using the boundary condition Pi(—1) = 1 and the continuity of Pi(t) across t = r — 1 we get 

+ + + forT<r-l, 

{ln(— t) — ln2-|-ln(r-|-2)} for T > r — 1. 

Using the symmetry between 1 and 2 we also get the same expression for P 2 {'^)- In terms of 
the physical time t we have 


InPi(T) = 


(3.13) 


PiW=/>2(i) = 


e (r-|-1-t-e *) ^'^(r-|-2)^^ for t < — ln(l — r), 
^(r-|-2)/2j for t > — ln(l — r). 


(3.14) 


Therefore, the life expectancy of Ci is 


U = 


t Pi{t) dt = / Pi{t)dt 


= {r + 2f^ B 


2 + r 


■,2K,0] -B 


1 — r 


2K,0 + 


2^^+^K 


where B{x]p,q) is the incomplete beta function, defined as 


B{x]p,q) = / P ^(1 —^dt = 


rx/{l-x) yp-i 


(3.15) 


(3.16) 


' y ' ' y (!+!/)« ■” 

the two expressions being related by the transformation t = y/{y + 1). In terms of the life 
expectancy T = l/2iF of a single isolated object, we have 


fi = (r + 2)^/'^ B f 


1 1 \ ^ 

2 + i’T’^) 


l-T 1 \ 

~’f’V 


(3.17) 
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Figure 7: The past light-come of C 2 at T 2 and the past light cone of Ci between Ti and Ti 
for Ti < T 2 — r. 


C 2 also has the same life expectancy. (3.17) is somewhat larger than T, but that is simply a 
result of our initial assumption that both objects exist at t = 0. If both objects had started 
at the same space-time point and then got separated following some specihc trajectories, then 
there would have been a certain probability that one or both of them will decay during the 
process of separation; this possibility has been ignored here leading to the apparent increase 
in the life expectancy. However, for realistic values of r and T, which corresponds to r << 1 
and T ~ 1, the ratio ti/T remains close to unity. 

Let us now turn to the computation of the joint survival probability T’i 2 (t^i, t^ 2 )- In this 
case the probability that the hrst object undergoes vacuum decay between Ti and Ti -|-5ti and 
the second object survives at least till T 2 is given by —(5 ti (cIPi 2 (t^i,' t 2 )/clTi). On the other 
hand the same probability is given by 77 x Pi 2 {xi, T 2 ) times the volume of the past light-come 
of Cl between Ti and Ti -|- 5ti, excluding the region inside the past light cone of C 2 at T 2 . The 
relevant geometry has been shown in Figs. 7, 8 and 9 for different ranges of Ti and T 2 . The 
results are as follows: 


1. For Ti < T 2 — r the geometry is shown in Fig. 7. In this case Ci at Ti (and hence the 
whole of the past light cone of Ci between Ti and Ti -|- (5ti) is inside the past light cone 
of O 2 at T 2 . Therefore, the decay probability is zero and we have the equation: 


ainPi2(Ti,T2) 

dxi 


0 for Ti < T 2 — r . 


(3.18) 


2. For T 2 — r < Ti < T 2 -|- r the geometry is as shown in Fig. 8. In this case Ci at Ti and C 2 
at T 2 are space-like separated. The volume of the past light cone of Ci between Ti and 
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T 



Figure 8: The past light-come of C 2 at T 2 and the past light cone of Ci between Ti and Ti 
for T 2 — r < Ti < T 2 -|- r. 


Ti -|- (5 ti outside the past light cone of C 2 at T 2 is given by 



PTI 

|(Tl+T2-r) 


da 


2 2 
Ti r - Ti - T 2 ■ 


This gives 


ainPi2(Ti,T2) 

d^i 



for T 2 — r < Ti < T 2 -|- r . 


(3.19) 


(3.20) 


3. For 0 < T 2 -|- r < Ti, the geometry is shown in Fig. 9. In this case C 2 at T 2 is inside the 
past light cone of Ci at Ti and there is no intersection between the past light cone of Ci 
between Ti and Ti -|- (5ti and the past light cone of C 2 at T 2 . Therefore, the volume of 
the past light cone of Ci between Ti and Ti -|- 5ti is given by 


and we have 


2 





ainPi2(Ti,T2) 

dxi 


2K — for T 2 -t- r < Ti < 0 . 


(3.21) 


(3.22) 


We can now determine Pi 2 (t^i,T 2 ) by integrating (3.18), (3.20), (3.22) subject to the boundary 
condition given in (3.9) 

-Pi2(~l5 ^^ 2 ) =-P2(t^2) =-Pi(t^2) 1 (3.23) 
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Figure 9: The past light cone of C 2 at T 2 and the past light cone of Ci between Ti and Ti + (5ti 
for T 2 + r < Ti < 0. 


and using the fact that -Pi 2 (t^i, ^^ 2 ) must be continuous across the subspaces dehned by Ti = 
T 2 ± r. The result of the integration is 


{ 2K {ln(—T 2 ) + ln(r + 2) — In 2} for Ti < T 2 — r , 

2iF{ln(-T2) + In(-Ti) - ln(r - Ti - T 2 ) + ln(r + 2)} , for T 2 - r < Ti < T 2 + r 
2iF{ln(—Ti) + ln(r + 2) — In 2} for T 2 + r < Ti < 0 . 

(3.24) 

Note that the result is symmetric under the exchange of Ti and T 2 even though at the inter¬ 
mediate stages of the analysis this symmetry was not manifest. 

Expressed in terms of physical time the above solution takes the form: 


I {(r -|- 2)/2}^^e“^^*2 for ti < — ln(r -1- , 

-^ 12 (^ 1 , ^ 2 ) = I (?" + -1- -f for — ln(r e~*^) < ti < — ln(e“*2 — r), 

[{(r -|- 2)/2}^^e“^^*^ for ti > — ln(e“*^ — r) . 

(3.25) 

If — r is negative then the third case is not relevant and in the second case there will be no 
upper bound on ti. Physically this can be understood by noting that in this case T 2 > —r and 
C 2 will never come inside the past light cone of Ci even when Ti reaches its maximum value 0. 

Our interest lies in computing the probability that at least one of the two objects survives 
till time t. Let us denote this by Puit). This is given by the sum of the probability that Ci 
survives till time t and the probability that C 2 survives till time t, but we have to subtract 
from it the probability that both Ci and C 2 survive till time t since this will be counted twice 
otherwise. This can be seen from the Venn diagram of two objects shown in Fig. 10. Therefore, 
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P{A UB) = P{A) + P{B) - P{A n B) 
Figure 10: Probability rule for N=2 using Venn Diagram. 


we have 

-Pi2(t) = -Pi(t) + -P2(t) — Puit, t) ■ (3.26) 


From this we can compute the probability that the last one to survive decays between t and 
t + 6t a.s 

-Stj^Pi2{t). (3.27) 

Therefore, the life expectancy of the combined system is given by 

d 

dtt—Pi2{t,t) = J dt {Pi{t) + P2{t) — Pi2{t,t)} , (3.28) 


where in the second step we have integrated by parts and used (3.26). Each of the hrst two 
integrals gives the result ii computed in (3.17). For the last integral since we have to evaluate 
Pi 2 iti,t 2 ) at ti = t 2 = t only the middle expression in (3.25) is relevant, and we get 


P^2{t,t)dt = I {r + 2f^e-^^\r + 2e-Y‘'^dt 


_t\-2K, 


= ir + 2Y^B 4K, -2K . 

2 + r 


(3.29) 
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Combining this with the result for ti given in (3.17) and replacing K by 1/2T we get 


tl2 


2{r + 2f/^ 


B 


1 1 
2 + r’T’° 


B 




2 + r’T’ 



+ T 


(1 — 


(3.30) 


We can now check various limits. First of all we can study the r —)■ 0 limit using the result 

5 ( 0 :; 2a, —a) ~ — (1 — x)“" , (3.31) 

a 

for X close to 1. This gives hmj._,.o t^i 2 = T. This is in agreement with the fact that if the two 
objects remain at the same point then their combined life expectancy is the same as that of 
individual objects. 

If on the other hand we take the limit of large T then, using the result 

B{x] a, /3) ~ — (3.32) 

a 

for small a, we get ti 2 — 3T/2. Therefore, the life expectancy of the two objects together is 3/2 
times that of an isolated object. This is consistent with the fact that if the inverse decay rate 
of individual objects is large then typically there will be enough time for the two objects to 
go out of each other’s horizon before they decay. Therefore, we can treat them as independent 
objects and recover the result (2.3). Mathematically this can be seen from the fact that when 
T is large and t ~ T then Pi 2 {t,t) given in the middle expression of (3.25) approaches 
which in turn is approximately equal to the square of Pi{t) given in (3.14). 


4 Vacuum decay in 3+1 dimensional de Sitter space 

In this section we shall repeat the analysis of §3 for 3+1 dimensional de Sitter space-time. 
Since the logical steps remain identical, we shall point out the essential differences arising in 
the two cases and then describe the results. 

The metric of the 3+1 dimensional de Sitter space is given by 

ds^ = —dt^ + e^^{dx‘^ + dy"^ + dz^) = r~'^(—dr^ + dx^ + dy^ + dz^), t = —. (4.1) 

There are of course various other coordinate systems in which we can describe the de Sitter 
metric, but the coordinate system used in (4.1) is specially suited for describing out universe. 
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with [x, y, z) labelling comoving coordinates and t denoting the cosmic time in which the 
constant t slices have uniform microwave background temperature. This form of the metric 
uses the observed flatness of the universe. The actual metric at present is deformed due to the 
presence of matter density, and also there is a lower cut-off on t since our universe has a hnite 
age of the order of the inverse Hubble constant. But both these effects will become irrelevant 
within a few Hubble time and we ignore them. In §5.3 we shall study these effects, but at 
present our goal is to get an analytic result under these simplifying assumptions. 


4.1 Isolated comoving object 

First consider the case of an isolated object. The calculation proceeds as in §3.1. However, in 
computing the volume of the past light cone in Fig. 3 we have to take into account the fact 
that for each a, the light cone is a sphere of radius (t — a). Since the coordinate radius of the 
sphere is (t — a) and the space-time volume element scales as 1/cr^ we get the volume of the 
past light cone of the object between t and t -|- (5t to be 


(5t 


— 47r(T— a) =-TTT OT. 

— rvn ^ 3 


(4.2) 


This replaces the right hand side of (3.4). Therefore, (3.5) takes the form 


with the solution 

Po{x) , 

(4.3) 


InFo(T^) = -7riFln(-T), 

(4.4) 


Po{t) = exp . 

(4.5) 

From this we can 

calculate the life expectancy of the isolated object to be 

/•oo 3 



^ = i = 4.A-■ 

(4.6) 


4.2 A pair of comoving objects 

The additional complication in the case of two objects comes from having to evaluate the 
contribution of the past light come of the hrst object between Ti and Ti -|- 5ti in situations 
depicted in Figs. 5 and 8. Let us consider Fig. 8 since Fig. 5 can be considered as a special 
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case of Fig. 8 with T 2 = —1. Now in Fig. 8 which occurs for T 2 — r < Ti < T 2 + r, the past 
light cone of Ci between Ti and Ti + lies partly inside the past light cone of C 2 . We need 
to subtract this contribution from the total volume of the past light cone of Ci between Ti 
and Ti + 5 ti, since the assumption that C 2 survives till T 2 rules out the formation of a bubble 
inside the past light cone of C 2 - Our goal will be to calculate this volume. 

Examining Fig. 8 we see that the intersection of the past light cones of Ci at Ti and C 2 at 
T 2 occur at T = a for cr < (ti + T 2 — r)/2. At a value of a satisfying this constraint, the past 
light cone of Ci at Ti is a sphere of coordinate radius ri = (ti — a) and the past light cone 
of C 2 at T 2 is a sphere of coordinate radius r 2 = (t 2 — a). The centers of these spheres, lying 
at the comoving coordinates of the two objects have a coordinate separation of r. A simple 
geometric analysis shows that the coordinate area of the part of the first sphere that is inside 
the second sphere is given by 

TT^ - (ri - rY} = —^(t 2 - Ti + r)(Ti + T 2 - r - 2a) . (4.7) 

Taking into account the fact that physical volumes are given by 1/cr^ times the coordinate 
volume we get the following expression for the volume of the past light cone of Ci between Ti 
and Ti + 5 ti that is inside the past light cone of C' 2 : 


„ ^(Tl+T2-r)/2 , 

-(t 2 - Ti + r) (5 ti / - cr)(Ti + T 2 - r - 2cr) 

r J-00 

271 (3r-Ti-3T2) 

— T2 - Ti + r 5 ti --— . 

3 r (r — Ti — T2)^ 


(4.8) 


As already mentioned the excluded volume in case of Fig. 5 can be found by setting Ti = t 
and T 2 = —1 in (4.8). 

We are now ready to generalize all the results of §3. Let us begin with (3.12). Its general¬ 
ization to the 3-1-1 dimensional case takes the form 

27iK 




3 

dvr a: 


2 (—1-1-r — t)( 3-I-3r — t) 

T r(T — r — 1)2 


if T < r — 1 
if T > r — 1 


(4.9) 


Its solution is given by 

^ In(-T) + g- ln(-T + r + 1) + + ln(r + 2) + 


InPi(T) = 


5r+6 

2r(r+2) 


(ln(— t) -|- ln(r -1- 2) — In 2 — 


2(r+2) 


if T > r — 1. 


if T < r — 1, 

(4.10) 
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ratio 



Figure 11: The figure showing the ratio ti/T for r = .0003, .001, .003, .01, .03, .1 and .3. For 
r < .003 the ratio is not distinguishable from 1 in this scale. 


Expressing this in terms of t using x = —e * and T = 3/(dvriF) we get 


n(() 



\ — — / \ — 

+ r + 1) r(^r + 2 )rexp 




t 1 r e-* 2(r + 1) 5r + 6 1 

T T \ 2r r(e“* + r + 1) 2r(r + 2) J 
for t < — ln(l — r), 


r + 2 


exp 


T 2T(r + 2)_ 


for t > — ln(l — r) 


(4.11) 

The same expression holds for the survival probability P 2 {t) of C 2 - From this we can find the 
life expectancy of Ci 


U = 


Pi{t)dt . 


(4.12) 


^'O 

As in the 1+1 dimensional case, p is slightly larger than T but this is simply due to the choice 
of initial condition that both observers are assumed to exist at f = 0. In Fig. 11 we have 
plotted the ratio P/T as a function of T for various values of r, and as we can see the result 
remains close to 1. More discussion on ii can be found below (4.23). 

Next we consider the generalization of (3.18)-(3.22). The analysis is straightforward and 
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we get the results 


ainPi2(Tl,T2) 

= 0 

d^i 

ainPi2(Ti,T2) 

2nK 

<9ti 

3 

ainPi2(Ti,T2) 

47rP 

9Ti 

3ti 


^ (r - Ti + T 2 )( 3 r - Ti - 3T2) 
Ti r (r - Ti - T2)2 

for T2 + r < Ti < 0 . 


for T2 — r < Ti < T2 + r 

(4.13) 


The solution to these equations, subject to the boundary condition Pi 2 (ti = —1, T 2 ) = P 2 (t^ 2 ) = 
Fi(t 2 ) is given by 


lnPi2(Ti,T2) = < 


4:TtK 

3 

AivK 

3 


3 


In 2 


if Ti < T2 — r 

T 1 +T 2 


ln(-T 2 ) + ln(r + 2) - - 

In(-Ti) + ln(-T 2 ) - ln(-Ti - T 2 + r) + - 

+ ln(r+ 2) + ^ 

In(-Ti) + ln(r + 2) - 2 ^ - In 2 


r(Ti+T 2 —r) 

if T2 - r < Ti < T2 + r 
if T2 + r < Ti < 0 


(4.14) 


In terms of the physical time, and T = 3/(47rP), this becomes 


{(r + 2)/2}Vr 


exp 


^2 


2 T(r + 2 ) 


if ti < — ln(r + e 


(r + 2 )^/'^(e + e 


Pi2(P, ^ 2 ) — < 


X exp 


p if 

1 


-*2 


+ r) T 


_T(r + 2 ) 2-*^^1+^2) 


1 . ^ f ^ 2 

(e-*i + e-*2) + 


2 Tr 


if 


Tr -g T-gO+fe 

ln(r + e“*^) < fi < — ln(e “*2 — r 


{(r + 2)/2}Vr 


exp 


2 T(r + 2 ) 


if fi > — ln(e — r) 


This gives 

Pi 2 (t, f) = (r + (^2e"* + r)~^ exp 


■—t — e * + 
T Tr 


T r 2 e* + re 


2t 


(4.15) 


(4.16) 


In terms of this, and the functions Pi = P 2 given in (4.11), we can calculate the probability 
P 12 of at least one of the two objects surviving till time t using 


Pi2it) — Pi(t) + P2it) — Puit, t) 


(4.17) 
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and the combined life expectancy of two objects using the analog of (3.28) 


ti2 — / Pi2{i)dt — / dt {Pi{t)P2{t) — — “iti — / dtPi2{t,t)- (4.18) 


For the integral of Pi 2 {t,t) one can write down an expression in terms of special functions 
as follows. Dehning y via 


, 2 + r 

2 + re* = - 

y 


(4.19) 


for r 7 ^ 0, we get 


dt Pi 2 {t,t) = [r{r + 2) [ dy y M — 


Now, using the result 


2y 

2 + r 


-l-l/T 


exp 


-•l y“-l(l - 

dy - = B{a,c — a)<hi(a, b, c; u, v) 


(2 + r)T 

(4.20) 


(4.21) 


(1 -uyf 

with Rec> Rea > 0, |m|< 1, B the beta function and <Fi the confluent hypergeometric series 
of two variables (Humbert series), we get 

1 l/T 


T 

dt Pi2{t,t) = — 


(r + 2) 


2 12 2 1 
QT(r+2) <J) ( _ X -I-1-1-• -- 

1 ' T, T’ T’2 + r’ (2 + r)T 


(4.22) 


$1 has a power series expansion 


ilia,b,c-,u,v)^ M<1 




(4.23) 


where (a)m = a(a + 1) ■ ■ ■ (a + m — 1). 

Unfortunately we have not been able to hud an expression for R = dt Pi{t) in terms of 
special functions. However, we can write down a series expansion for this that will be suitable 
for studying its behaviour for small r. The integral of (4.11) from t = — ln(l — r) to oo is 
straightforward and yields 


T (1 - r)^/^(r + 2)^/^2-^/^exp 


r 

2T(r + 2) 


(4.24) 
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The integral of (4.11) from f = 0 to — ln(l — r) can be analyzed by making a change of variable 
from t to y via e~^ = {1 — yr). In terms of this variable the integral can be expressed as 


1 / \ 




exp 


2y‘^r{r + 1) f 


T(2 + r)3 V 2 + r 


-1 


exp 


Using series expansion of the second and third terms in the integrand we get 


2 ~\ 
yr 

2r(2 + r)2_ 

(4.25) 


E 

r72,n=0 


m\ n\ 


1\ 2V™+^+i(r + l)' 

’ T-(2 + r)3- 


dyy 


,m+2n 


yr 

2 + r 


-n-ljT 


exp 


yr 


_2T(2 + r)2_ ■ 


(4.26) 


The integral over y can be expressed in terms of <l>i using (4.21). Adding (4.24) to this we get 


ti = T(1-r)^/'^(r+ 2)^/^2-i/^exp 
1 1 


m! n! m + 2n + 1 I T J ^ ^ 

r72,n=0 / m 


2T(r + 2)_ 

r2nj.m+n+l^^ _j_ 


Tn{2 + r^n 


4.. I m + 2n + 1. n + 1. ,n + 2n + 2; 


(4.27) 


It can be checked using (4.11), (4.12), (4.16) and (4.18) that for r —0 we get U 2 A 1 — 1 
and for T —)■ cx) we get fi 2 /ti = 3/2. The values of ‘gain’ = fi 2 /ti for different values of r have 
been plotted against T in Fig. 1. 


4.3 The case of small initial separation 

Since from practical considerations the small r region is of interest, it is also useful to consider 
the expansion of ii 2 /ii for small r. For this we have to analyze the behaviour of + as well as 
that of dt P 12 (t,t) for small r. Let us begin with ti given in (4.27). It can be easily seen 
that this is given by T + 0{r) with the contribution T coming from the hrst term. However, 
the contribution from J^dtPi 2 (t,t) has a more complicated behaviour at small r. This is 
related to the fact that in the r —)• 0 limit the fourth argument of <Fi in (4.22) approaches 1, 
and in this limit the series expansion (4.23) diverges. To study the small r behaviour we shall 
go back to the original expression for Puit, t) given in (4.16). We change variable to n = e~^ jr 
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and write 


I t'^jv J 

dtPi2{t,t) = (r + 2)i/V/^("+2Vi/^ / —(2n + l)-'/V/'^ exp 

./n V 


= (r + 2)VV/^(^+2)^i/^ n 


2/T 


{ 2 v + 1) 


-l/T 


T{2v + 1) 

V 


-( 2 v 


,-l/T 


exp 


1 


+ (r+ 2)i/^2-^/'^ 


exp 


exp 

1 


r(2n + l)_ 

1 


_T(r + 2) 2T 


T, 

(4.28) 


where in the last step we have snbtracted an integral from the original integral and compen¬ 
sated for it by adding the explicit resnlt for the integral. This snbtraction makes the integral 
convergent even when we replace the npper limit 1/r by oo. Taking the small r limit we get 




(2t. + 1) 


-l/T 


exp 


V 1 

'T( 2 v + 1 ) ^ ^ 


- (2n)-^/^ 
(4.29) 


+ T + 0(r). 

Combining this with the earlier resnlt that for r —)■ 0, fi ~ T -|- 0(r) and using (4.18) we get 

(4.30) 


j. 1 /•oo 

-f=2-^ dtPi2{t,t) = l + A{T)r^/^, 
Jo 


where 

poo 

A{T)=T-^ 2 ^/'^ dvv-'+^^^ 


(2v) 


-l/T _ 


{ 2 v + 1) 


-l/T 


exp 


V 1 

' t {2 v + 1) ^ if 


(4.31) 


The numerical values of A{T) are moderate ~ for example ^(5) ~ 0.439 and 44(10) ~ 0.457. 
A plot of A{T) as a function of T has been shown in Fig. 12. The l/T exponent of r shows 
that even if we begin with small r, for moderately large T (say T ~ 5) we can get moderate 
enhancement in life expectancy. 


5 Generalizations 

In this section we shall discuss various possible generalizations of our results. 

5.1 Multiple objects in de Sitter space 

We shall begin by discussing the case of three objects Ci, C 2 and C 3 placed at certain points 
in 3 -I- 1 dimensional de Sitter space-time and analyze the probability that at least one of them 
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Figure 12: The coefficient of the term in the expression for tyijtx as a function of T. 

will survive till time t. Let -Pi 23 (fi,^ 2 ,^ 3 ) denote the probability that Ci survives till time ti, 
C 2 survives till time t 2 and C 3 survives till time t^. Similarly Pij{ti,tj) for 1 < < 3 will 

denote the probability that Ci survives till ti and Cj survives till tj and Pi{ti) will denote the 
probability that Ci survives till ti. All probabilities are dehned under the prior assumption 
that all objects are alive at t = 0. These probabilities can be calculated by generalizing the 
procedure described in §3 and §4 by constructing ordinary differential equations in one of the 
arguments at hxed values of the other arguments. The geometry of course now becomes more 
involved due to the fact that the past light cone of one object will typically intersect the past 
light cones of the other objects which themselves may have overlaps, and one has to carefully 
subtract the correct volume. But the analysis is straightforward. 

The quantity of direct interest is the probability P^sit) that at least one of the objects 
survives till time t. With the help of the Venn diagram given in Fig. 13 we get 

Pl23{t) = + P2{t) + Psii) — Pl2{t, t) — Pl3{t, t) — P23{t, t) + Pl23{t, t, . (5.1) 

Using this we can calculate the life expectancy of the combined system as 

d ~ ~ 

— J dtt—Pi23it) = J dtPi23{t)- (5-2) 
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Figure 13: The Venn diagram illustrating that the survival probability of one of A, B or C, 
denoted by P{A U S U (F), is given by P{A) + P{B) + P(C') - P{A fl P) - P{B fl C) - P{A n 
C) + P{AnBf\C). 

The generalization to the case of N objects is now obvious. The relevant formula is 

N N N 

Pi2-N{t) = [j^Pi{t) - '^Pij{tA) + ^ Pijk{tAA )^—• • • ,^)) (5.3) 

j=l i<j i<j<k 

where Pjj... 4 (tj^, • • are again computed by solving ordinary differential equations in one 
of the variables. Once Pi 2 . .Af(t) is computed we can get the life expectancy of the combined 
system by using 

tl2-N = / Pi2-N{t)dt . (5.4) 

Jo 

5.2 Realistic trajectories 

Another generalization involves considering a situation where multiple objects originate at 
the same space time point and then follow different trajectories, eventually settling down at 
different comoving coordinates. This has been illustrated in Fig. 14. This represents the 
realistic situation since by dehnition different civilizations of the same race must originate at 
some common source. We can now generalize our analysis to take into account the possibility 
of decay during the journey as well. Eqs.(5.3) and (5.4) still holds, but the computation of 
Pii-ikitii, ■ ■ ■ ti^) will now have to be done by taking into account the details of the trajectories 
of each object and the overlaps of their past light cones. The principle remains the same, and 
we can set up ordinary differential equations for each of these quantities. The only difference 
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Figure 14: Multiple objects originating from the same space-time point. Different dashed lines 
represent the trajectories followed by different objects. 

is that the spatial separation between the i-th object at t = t* and the j-th object at t = Xj 
will now depend on Tj and Xj according to the trajectories followed by them. 

This analysis can be easily generalized to the case where each of the descendant objects in 
turn produces its own descendants which settle away from the parent object and eventually go 
outside each other’s horizon due to the Hubble expansion. If this could be repeated at a rate 
faster than the vacuum decay rate then we can formally ensure that some of the objects will 
survive vacuum decay [16]. However, since within a few Hubble periods most of the universe 
will split up into gravitationally bound systems outside each other’s horizon, in practice this 
is going to be an increasingly difficult task. 


5.3 Matter effect 


A third generalization will involve relaxing the assumption that the universe has been de Sitter 
throughout its past history. While de Sitter metric will be a good approximation after a few 
Hubble periods, within the next few Hubble periods we shall still be sensitive to the fact that 
the universe had been matter dominated in the recent past and had a beginning. This will 
change the form of the metric (4.1) to 

ds^ = —dt^ -I- a{tf‘{dx^ -|- dy‘^ + dz^) (5.5) 


where a{t) is determined from the Friedman equation 


1 da 
a dt 



(5.6) 
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in the convention that the value of a is 1 today and pA and pm are the energy densities 
due to cosmological constant and matter today. Since we have chosen the unit of time so 
that the Hubble parameter in the cosmological constant dominated universe is 1, we have 
^/SttGpa/S = 1. Dehning® 

c = Pm/P a - 0.45 

we can express (5.6) as 


1 da r~ ^ 

—- = V1 + ca-3 . 
a dt 


Let T be the conformal time defined via 


dx = dt/a{t) 

with the boundary condition T —)■ 0 as f —)■ oo. Then (5.6) takes the form 

1 da n -^ 

^^ = vl + ca 
dr 


whose solution is 


T = 


db 






a/1 + 


dv 1/114 c 

^~~~a Is’ 2’ 3’“^ 


This implicitly determines a as a function of t. The metric is given by 

ds^ = a{r)‘^{—dr'^ + dx'^ + + dz^). 


(5.7) 

(5.8) 

(5.9) 

(5.10) 

(5.11) 

(5.12) 


Using the experimental value c ~ 0.45 we get that t —)■ Tq ~ —3.7 as a —)■ 0, showing that 
the big bang singularity is at t ~ —3.7.^ We also have that at a = 1, t ~ —0.95. This is 
not very different from the value T = — 1 for pure de Sitter space-time with which we have 
worked. However, we shall now show that the decay rate in the matter dominated epoch 
of the universe differs significantly from that in the cosmological constant dominated epoch. 
The decay rate of an isolated observer at some value of the conformal time t is given by the 
following generalization of (4.2), (4.3): 


d 

dr 


InPo(T^) 


-ArK 

-AnK 


da a(ar ("^ — a) 


"To 

/.a(T) 


db 


'b=o ft^Vl + cb ^ 



(I 14 
V3’2’3’ 



(5.13) 


®We use cosmological parameters given in [18]. 

^Of course close to the singularity the universe becomes radiation dominated but given the short span of 
radiation dominated era we ignore that effect for the current analysis. 
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where in the second step we have changed the integration variable from a to b = a{a). From 
this we can compute the decay rate: 


m 


p ^_ 

dt ° a{t) dr 
AttK db 

a(t) 7^=0 a /1 + cb-^ 


\nPo{r) 

b^ |T(t) + 


1 

b 


2F1 


1 1 4 
3’ 2’ 3’ 




(5.14) 


Using the information that today a = 1 and t ~ —0.95 we get 

,,, AnK 0.067 

11(f) Itoday^ — X 0.067 ~ . (5.15) 

This is lower than the corresponding rate T~^ in the de Sitter epoch by about a factor of 15. 
The growth of the decay rate with scale factor has been shown in Fig. 2. 

Our analysis of §4 can now be repeated for two or more observers and also for general tra¬ 
jectory discussed in §5.2 with this general form of the metric to get more accurate computation 
of the life expectancy. These corrections will be important if T ~ 1 and the decay takes place 
within a few Hubble period from now. On the other hand if T is large (say ~ 10) then the 
decay is likely to take place sufficiently far in the future by which time the effect of our matter 
dominated past will have insignificant effect on the results. 

The fact that the decay rate increases with time till it eventually settles down to a constant 
value in the de Sitter epoch has some important consequences: 


1. We have already seen from (5.15) that the decay rate today is about 15 times smaller 
than the decay rate in the de Sitter epoch. Eq.(5.14) for a{t) = 2 shows that even when 
the universe will be double its size compared to today, the decay rate will remain at about 
27% of the decay rate in the de Sitter epoch. Since most of the journeys to different parts 
of the universe ~ if they take place at all - are likely to happen during this epoch, we see 
that the probability of decay during the journey will be considerably less than that in the 
final de Sitter phase. This partially justifies our analysis in §4 where we neglected the 
probability of decay during the journey. This also shows that if we eventually carry out 
a detailed numerical analysis taking into account the effect discussed in §5.2, it should 
be done in conjunction with the analysis of this subsection taking into account the effect 
of matter. 


2. It is also possible to see from (5.14), (5.15) (or Fig. 2) that the decay rate in the past 
was even smaller than that of today. If D{t) denotes the decay rate at time t dehned in 









(5.14), then the average decay rate in onr past can be dehned as 


1 

- / D{t)dt, 

Jo 

where ti denotes the cnrrent age of the nniverse given by 

V dt 


(5.16) 


, f da 
ti= da { — 


da 


lo ay/T+ 


ca 


-3 


0.79. 


(5.17) 


In physical units ti is about 1.38 x lO^'^ years. The evaluation of (5.16) can be facilitated 
using the observation that D{t)dt is K times the volume enclosed between the past light 
cones of the object at times t and t + dt. Therefore, dtD{t) must be K times the total 
volume enclosed by the past light cone of the object at ti. This can be easily computed, 
yielding 




^ dh r 1 

T + cb-^ \ b 


2F1 
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(5.18) 


For c ~ 0.45 and t given by today’s value —0.95 this gives 

D(t)dt = . (5.19) 

Using (5.16), (5.17) we get the average decay rate to be 0.018/T. This is about 3.7 times 
smaller than the present decay rate given in (5.15) and 56 times smaller than the decay 
rate 1/T in the de Sitter epoch. 



Integrating the equation dP^/dt = —D(t)PQ(t) we get 


InFo(U) 



0.014 

T 


(5.20) 


Requiring this to be not much smaller than —1 (which is equivalent to requiring that the 
inverse of the average decay rate (5.16) be not much smaller than the age of the universe 
ti) gives T ~ 0.014. This is much lower than what one might have naively predicted by 
equating the lower bound on T to the age of the universe i.e. T ~ ~ 0.79. Recalling 

that the unit of time is set by the Hubble period in the de Sitter epoch which is about 
1.7 X 10^° years, the bound T ~ 0.014 translates to a lower bound of order 2.5 x 10® 
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years. Since the current decay rate is about 15 times smaller than that in the hnal de 
Sitter epoch, we see that the lower bound on the current inverse decay rate is of order 
3.7 X 10® years. This is comparable to the period over which the earth is expected to be 
destroyed due to the expanded size of the Sun. 

3. Finally we note that the above analysis was based on the assumption that the bubbles 
continue to nucleate and expand in the FRW metric at the same rate as they would do 
in the metastable vacuum. This will be expected as long as the matter and radiation 
density and temperature are small compared to the microscopic scales involved in the 
bubble nucleation process, e.g. the scale set by the negative cosmological constant of 
the vacuum in the interior of the bubble. Some discussion on the effect of cosmological 
space-time background on the bubble nucleation / evolution can be found in [19,20]. 

6 Discussion 

We have seen that the result for how much we can increase the life expectancy by spreading 
out in space depends on the parameters r and K, which in turn are determined by the Hubble 
parameter of the de Sitter space-time, the initial spread between different objects and the 
inverse decay rate. Therefore, the knowledge of these quantities is important for planning our 
future course of action if we are to adapt this strategy for increasing the life expectancy of the 
human race. In this section we shall discuss possible strategies for determining / manipulating 
these quantities. 

We begin with the Hubble expansion parameter H. This is determined by the cosmological 
constant which has been quite well measured by now. Assuming that the current expansion 
rate is of order 68Km/sec/Mpc and accounting for the fact that the cosmological constant 
accounts for about 69% of the total energy density we get ~ 1.7 x 10^® years. Future 
experiments will undoubtedly provide a more accurate determination of this number, but 
given the uncertainty in the other quantities, this will not signihcantly affect our future course 
of action. Of course we may discover that the dark energy responsible for the accelerated 
expansion of the universe comes from another source, in which case we have to reexamine the 
whole situation. 

Next we turn to the initial separation between different objects which determine the value 
of r. Since in order for the Hubble expansion to be effective in separating the objects they 
have to be unbound gravitationally, a minimum separation between the objects is necessary for 
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overcoming the attractive gravitational force of the home galaxy. For example the size of our 
local gravitationally bound group of galaxies is of order 5 million light years which correspond 
to r ~ 3 X 10““^. The question is whether larger values of r can be accessed. An interesting 
analysis by Heyl [21] concluded that by building a space-ship that can constantly accelerate / 
decelerate at a value equal to the acceleration due to gravity, we can reach values of r close 
to unity in less than 100 years viewed from the point of view of the space-traveller. Of course 
this will be close to about 10^° years viewed from earth, and roughly the reduction of time 
viewed from the space-ship can be attributed to the large time dilation at the peak speed of 
the space-ship reaching a value close to that of light. However, this large time dilation will 
also increase the effective temperature of the microwave background radiation in the forward 
direction and without a proper shield such a journey will be impossible to perform. If one 
allows a maximum time dilation of the order of 100 then the microwave temperature in the 
forward direction rises to about the room temperature. Even then we have to worry about 
the result of possible collisions with intergalactic dust and othe debris in space. Even if these 
problems are resolved, we shall need a time of order 10® years from the point of view of the 
space-ship to travel a distance of order 10^° light years. Even travelling the minimum required 
distance of order 10^ light years will take 10^ years in such a space-ship. Such a long journey 
in a space-ship does not seem very practical but may not be impossible. 

Another interesting suggestion for populating regions of space-time which will eventually be 
outside each other’s horizon has been made by Loeb [22] . Occasionally there are hypervelocity 
stars which escape our galaxy (and the cluster of galaxies which are gravitationally bound) 
and so if we could hnd a habitable planet in such a star we could take a free ride in that planet 
and escape our local gravitationally bound system. In general of course there is no guarantee 
that such a star will reach another cluster of galaxies where we could spread out and thrive, 
but some time we may be lucky. It has been further suggested in [23,24] that the merger of 
Andromeda and the Milky Way galaxies in the future [25] could generate a large number of 
such hypervelocity stars travelling at speeds comparable to that of light and they could travel 
up to distances of the order of 10® light years by the time they burn out. This could allow us 
to achieve values of r of order 10“^ or more. 

Let us now turn to the value of T or equivalently the decay rate of the de Sitter vacuum 
in which we currently live. This is probably the most important ingredient since we have 
seen that for T ~ 1 we do not gain much by spreading out, while for large enough T we can 
achieve the maximum possible gain, given by the harmonic numbers, by spreading out even 
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over modest distances of r ~ 10“^. At the same time if T is so large that it exceeds the period 
over which galaxies will die then vacnnm decay may not have a significant role in deciding onr 
end and we should focus on other issues. For this reason estimating the value of T seems to be 
of paramount importance. Unfortunately, due to the very nature of the vacuum decay process 
it is not possible to determine it by any sort of direct experiment since such an experiment will 
also destroy the observer. It may be possible in the future to device clever indirect experiments 
to probe vacuum instability without actually causing the transition to the stable vacuum, but 
no such scheme is known at present. 

At a crude level the current age of the universe - which is about 0.79 times the asymptotic 
Hubble period in the cosmological constant dominated epoch - together with the assumption 
that we have not been extremely lucky to survive this long, suggests that the inverse decay 
rate of the universe is ~ 0.8. However, (5.18) shows that this actually gives a lower bound of 
T ~ 0.014 reflecting the fact that the decay rate in the de Sitter epoch will be about 56 times 
faster than the average decay rate in the past. If we allow for the possibility that we might 
have been extremely lucky to have survived till today, then we have indirect arguments that 
lower the bound by a factor of 10 [13]. Clearly these rates are too fast and if T really happens 
to be less than 1, then there is not much we can do to prolong our collective life. 

Is there any hope of computing T theoretically? Unfortunately any bottom up approach 
based on the analysis of low energy effective field theory is insufficient for this problem. The 
reason for this is that the vacuum decay rate is a heavily ultraviolet sensitive quantity. Given 
a theory with a perfectly stable vacuum we can add to it a new heavy scalar field whose effect 
will be strongly suppressed at low energy, but which can have a potential that makes the 
vacuum metastable with arbitrarily large decay rate. For this reason the only way we could 
hope to estimate T is through the use of a top down approach in which we have a fundamental 
microscopic theory all of whose parameters are fixed by some fundamental principle, and then 
compute the vacuum decay rate using standard techniques. In the context of string theory this 
will require finding the vacuum in the landscape that describes our universe. Alternatively, in 
the multiverse scenario, we need to carry out a statistical analysis that establishes that the 
overwhelming majority of the vacua that resemble our vacua will have their decay rate lying 
within a narrow range. This can then be identified as the likely value of the decay rate. There 
have been attempts in this direction [26-29], but it is probably fair to say that we do not yet 
have a definite result based on which we can plan our future course of action. 
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A Decay rate for equation of state p = w p 

In this appendix we shall compute the growth of decay rate with time for a general equation 
of state of the form p = w p with w > —1. In this case the p and a are related as 




(A.l) 


for some constant po- As a result the dependence of a on f and t are determined by the 
equations 
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The solutions to these equations are 
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As a ranges from 0 to oo, both t and t also range from 0 to oo. 

We can now compute the decay rate using the hrst equation of (5.13): 

D{t) = --^lnPo('r) 
a dr 

= -AttK f daa{ay{x — ay, 

® Jo 


(A.5) 


where a(cr) denotes the scale factor at conformal time a. Changing integration variable to 
b = a(a), which due to (A.4) corresponds to 
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we can express (A.5) as 
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Using the relation between a and t given in (A.4), this can be expressed as 

{w + lf 
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